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Abstract 

We derive the structural relations between the Mellin transforms of weighted Nielsen inte- 
grals emerging in the calculation of massless or massive single-scale quantities in QED and 
QCD, such as anomalous dimensions and Wilson coefficients, and other hard scattering 
cross sections depending on a single scale. The set of all multiple harmonic sums up to 
weight five cover the sums needed in the calculation of the 3-loop anomalous dimensions. 
The relations extend the set resulting from the quasi-shuffle product between harmonic 
sums studied earlier. Unlike the shuffle relations, they depend on the value of the quanti- 
ties considered. Up to weight w = 5, 242 nested harmonic sums contribute. In the present 
physical applications it is sufficient to consider the sub-set of harmonic sums not containing 
an index i = — 1, which consists out of 69 sums. The algebraic relations reduce this set to 
30 sums. Due to the structural relations a final reduction of the number of harmonic sums 
to 15 basic functions is obtained. These functions can be represented in terms of factorial 
series, supplemented by harmonic sums which are algebraically reducible. Complete an- 
alytic representations are given for these 15 meromorphic functions in the complex plane 
deriving their asymptotic- and recursion relations. A general outline is presented on the 
way nested harmonic sums and multiple zeta values emerge in higher order calculations of 
zero- and single scale quantities. 



1 Introduction 

Finite harmonic sums [1-3] and associated to them, by a Mellin transform, Nielsen-type inte- 
grals [4], emerge in perturbative calculations of massless and massive single scale problems in 
Quantum Field-Theory. The Wilson coefficients and anomalous dimensions in deeply inelastic 
scattering to three loops [5], the massive quark Wilson coefficients in the limit Q 2 3> m 2 [6,7], as 
well as the Wilson coefficients for the Drell-Yan process, hadronic Higgs production in the heavy 
mass limit, the time-like coefficient functions for parton fragmentation into hadrons, cf. [8], the 
soft- and virtual corrections to Bhabha scattering [9] , and various processes more belong to this 
class. At 2-loop order the respective expressions in terms of harmonic sums were given in [10] 
using the algebraic and structural relations between these quantities of weight w < 4. Due to the 
complexity of higher order calculations the knowledge of as many as possible relations between 
the finite harmonic sums is of importance to simplify the calculations and to obtain compact 
analytic results. It turns out that, up to isomorphisms, the basic functions, which express the 
physical quantities mentioned above, are always the same and form a basis of an algebra, over 
which the observables can be expressed. These basic functions therefore form central objects in 
the description of a wide range of physical quantities. 

To derive the basic functions we first consider the multiple finite harmonic sums S ait ... ,a„(N). 
They provide a unique language to express single scale quantities, which emerge in field theoretic 
calculations. They are defined by 



Here, are positive or negative integers and N is a positive even or odd number. One calls n 
the depth and YX=\ |ofc| the weight of a harmonic sum. Harmonic sums are associated to Mellin 
transforms of real functions f(x) 



Here M denotes the Mellin-transform. The functions f ai ,...,a n {x) are usually linear combinations 
of harmonic polylogarithms Hi, Ummm ^(x) [11-13], weighted by the factor 1/(1 ±x). f(x) may be 
distribution-valued [14]. One example for this class is 



The number of harmonic sums grows exponentially with the weight w and amounts to 3™- 1 
for all sums of weight w' < w. The finite harmonic sums are not independent objects but obey 
relations. In view of their rapidly growing number for higher weight it is worthwhile to reveal all 
their relations to express all harmonic sums in terms of a small set of basic elements only. In a 
foregoing paper [15] we investigated the shuffle relations [16] between general finite alternating 
or non-alternating harmonic sums. 1 There all shuffle relations for harmonic sums up to depth 
and weight w = 6 were derived in general form and the counting relations for the number of 
basis elements spanning the algebra of harmonic sums after the shuffle relations are used have 

1 Various mathematical terms used in the following were defined in Ref. [15] to which we refer readers from 
the physics community who are not yet used to this terminology. 
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been given. For a certain class of indices {a±, . . . , a*,}, where rij\j=i,,, q elements are equal and 
n = J2 9 j=i n ji the total number of harmonic sums is 



n(ni, ... ,n q ) 



ri] \ . . . n q \ 



(1.4) 



The shuffle relations allow to express these sums in terms of 



Z n (ni,... ,n 9 ) = — A*(d) 

d\iii 




(1.5) 



(f)'---(^)' 



basis elements and polynomials out of objects of lower weight. Here Eq. (1.5) is the 2nd Witt 
formula, cf. [17], and fi(k) denotes the Mobius function [18]. l n also counts the number of 
Lyndon words [19] associated to {ai, . . . , a^}, since all harmonic sums belonging to this index 
set are freely generated by the Lyndon words as a consequence of Radford's theorem [20]. 

In the limit iV — > oo harmonic sums turn into multiple zeta values (MZVs) [21] 2 . Since 
multiple zeta values obey many more relations, see e.g. [22,23], than the shuffle relations the 
question arises which other relations may hold for multiple harmonic sums. 3 A series of relations 
emerges due to the integral representations of the harmonic sums which were derived in Ref. [2]. 
Furthermore, some sums have a completely symmetric index set and do therefore decompose into 
polynomials of single harmonic sums. This also implies relations between the Mellin transforms 
of certain Nielsen-type integrals weighted with the denominators 1/(1 ± x). 

In the present paper we will not limit the consideration to harmonic sums with iV G N but 
will allow for N <E C and derive the complex analysis for these quantities. Along with this, 
new relations between these objects will emerge. In this generalized case we consider rational 
fractions of N, N/q, and allow to differentiate for the index N, which both yield new relations. 
The harmonic sums turn out to be meromorphic functions. They can be represented through 
factorial series [27,28] up to polynomials of Si(N) and harmonic sums of lower degree. The 
analytic continuation of S\(N) is given by Euler's di-gamma function ip(z). 

The weight w = 1 harmonic sums S±k(N) may all be traced back to the ^—function for 
iV G C. At the beginning of the 20th century Nielsen [27] studied the respective functions 
belonging to the class of weight w = 2, from which the corresponding analytic continuations can 
be derived. Due to the fact, that harmonic sums, with the exception of factors of S±(N), can 
be represented as factorial series, allows to derive their analytic continuations. The recursion 
relations for (N + 1) — > N, N e C, and the asymptotic representation for \N\ — > oo are 
obtained in analytic form. The singularities are poles located at the non-positive integers. 
Due to the presence of factors Si(N) in some cases the respective sums behave oc \n k (N) as 
\N\ — > oo. 4 Since we derive the analytic continuations in analytic form, also all derivatives are 

2 If oi = 1, the corresponding multiple sums diverge for TV — > oo but can be considered symbolically. Actually 
they are expressible in terms of a polynomial of the (divergent) harmonic series oo = X)fc°=i( 1 /^) an( i the finite 
basis elements of multiple zeta values, see [3,15,24]. 

3 The relations between the alternating MZVs to weight w = 12 and non-alternating MZVs to weight w = 24 
based on the shuffle-, stuffle-, and new relations have been determined in [24] recently. Generalizations of finite 
harmonic sums and harmonic polylogarithms have been given for other applications, cf. [25] and [26]. 

4 Precise analytic continuations of the basic functions [2, 29] based on semi-numerical representations were 
given in [30,31]. Here we made use of the MINIMAX method [32,33]. This method has also been applied to derive 
the analytic continuation for the heavy flavor Wilson coefficients up to 2-loop order [34]. For another proposal 
for the analytic continuation of harmonic sums to N <G R, for which some simple examples were presented, 
cf. [35] . For other effective parameterizations see [36] . 
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easily determined. In the definition of the basic functions we chose representations in the past, 
in which the functions to be Mellin-transformed may posses a branch point at x — 1. To express 
the Mellin transforms M[/(x)](iV) in terms of factorial series it is required that the functions 
f(x) are analytic at x — 1, which can be obtained using well-known mirror relations for the 
corresponding Nielsen integrals or harmonic polylogarithms transforming x <-> (1 — x). 

The paper is organized as follows. In Section 2 a brief general outline is presented on the way 
finite nested harmonic sums emerge in higher order calculations for single scale processes. Some 
basic definitions are given in Section 3. Classes of harmonic sums are discussed in Section 4. 
In sections 5 to 9 we consider the Mellin transforms being associated to the sums of weight 
w = 1 to w = 5 in detail and derive the structural relations for the harmonic sums. 5 Here 
we will first refer to the set of Mellin transforms chosen formerly [2,15,30]. The basis found 
in this way will be reordered by a linear transformation into another one which is consistent 
with the representation in terms of factorial series. Section 10 contains the conclusions. The 
algebraic relations to weight w = 5 used are given in Appendix A. Appendix B contains the 
relations necessary to express the basic functions in the complex plane. Some useful integrals 
which emerge in this context are given in Appendix C. 

2 The Emergence of Harmonic Sums 

In the following we give a brief outline on the way finite nested harmonic sums emerge in higher 
order calculations for single scale processes. These processes contain a scale, such as a momentum 
fraction z, which relates a parton momentum p = zP collinearly to the hadron momentum P 
in the massless limit, < z < 1. Examples for this case are space- and time-like splitting and 
coefficient functions, cf. [5,8]. The representation also applies to massive calculations in the limit 
Q 2 /m 2 3> 1, cf. e.g. [6,7]. Likewise this variable may be viewed as the ratio of two Lorentz- 
invariants z = t/s, as the case for the soft and virtual corrections to Bhabha-scattering [9], or 
the initial state radiation to e + e~ annihilation [38], where the respective ratio is z = s' / s with s' 
the cms energy of the virtual gauge boson exchanged. In the latter two processes the electrons 
are massive but are dealt with in the limit s/m 2 e ^> 1. s, s' and t are Mandelstam- variables, i.e. 
Lorentz-products of linear combinations of 4- momenta. The factions z do always obey < z < 1. 

There are two principle approaches, which can be followed in these calculations. Either one 
uses the light-cone expansion [39] for space-like and the cut-vertex method [40] in the time-like 
case or one follows a partonic description in the massless case, resp. the case Q 2 /m 2 ^> 1, with 
Q 2 the relevant bosonic virtuality exchanged. The light-cone expansion and cut-vertex method 
lead to Mellin moments immediately. Here, furthermore, the corresponding current crossing 
relations have to be observed, which single out even or odd integer moments, cf. [41], depending 
on the physical process. In the partonic approach one may form these moments accordingly, 
transforming to Mellin space since the variables t/s, resp. s' / s G [0, 1] in the limit mentioned 
above. Further on we will consider the representations in Mellin space only. 

The no-scale quantities can be obtained as a special case out of single scale quantities setting 
N = const, or N — > oo, respectively. If in the former case, the representation appears in terms 
of nested harmonic sums, partly decorated with multiple (-values, the latter are just multiple 
(■-values. Here we always refer to the alternating harmonic sums and multiple (-values, as the 
present level of physics calculations for single scale processes contains these quantities only. These 
structures are obtained by iterated integrals over the alphabet 1/x, 1/(1 — x), 1/(1 + x) and 

5 The structural relations of w = 6 are presented in Ref. [37]. 
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their Mellin transforms. Going to even higher orders, or including a scale more, this alphabet will 
be extended, cf . [42] , which leads to similar structures. The role of the harmonic polylogarithms 
[12] will be taken by the polylogarithms over the extended alphabet. The multiple (^-values are 
then generalized by the values of the new polylogarithms at x — 1, resp. the symbolic polynomial 
in case the divergent symbol <jq appears. The nested harmonic sums are generalized into linear 
combinations of Mellin transforms of the new polylogarithms, which one may call generalized 
harmonic sums. 

Let us now describe how these structures emerge. We consider a renormalizable Quantum 
Field Theory in D = n + e space-time dimensions. In case of the Standard Model and its 
parts, Quantum Electrodynamics and Quantum Chromodynamics, the dimension is n = 4. The 
propagators and vertices of the Feynman diagrams are combined with the help of Feynman 
parameters, 

1 T(a + (3) f 1 dx x a ~ l (1 -xf- 



f+ dx [L-xr^ 

Jo \xA + (1 - x)B]<*+f> ' 1 ' ' 



A a ■ BP r(a)r(/3) Jo [xA + (1 - x)B] a +< 3 

This allows to re-express the otherwise difficult angular integrals. In this way the loop D- 
momenta are synchronized. The Wick-rotation is performed, which leads to Euclidean integrals 
for the loop momenta in D = n + e dimensions. The numerator integrals may contain operator 
insertions like 

(A.ki) m , (2.2) 

and several finite sums over these terms in m, which are bounded by the Mellin variable N. The 
operator insertion (2.2) in momentum space stem from local composite operators emerging in 
operator product- or light cone expansions. Here hi denotes a loop momentum, m an integer 
and A a light-like .D-vector with A. A = 0. The dot denotes the Minkowski-product, which turns 
into the scalar product through Wick-rotation. The trivial angular integrals and radial integrals 
over the loop momenta can now be performed, which yield a factor 

r(rai + r\E) . . . T(n k + r k e) 



T(m 1 + qie) . . .r(mj + q t e) 



(2.3) 

ii,m,j e Z,ri,qj £ Q 



leaving the non-trivial Feynman parameter integrals. In carrying out theses integrals scalar and 
tensor integrals have to be performed. 

The Feynman parameter integrals contain the Mellin variable in the numerator only. The 
denominator is given by the Kirchhoff-polynomial of the respective graph consisting out of k 
monomials raised to a real non-integer power. A possible next step consists in expressing the 
Feynman-parameter integrals in terms of Mellin-Barnes integrals [43] for each of these monomials, 

- =— da A" B~ q ~ a a >\y q ^ a > . (2.4) 

(A + B) q 2TxiJ 1 -ioo T(q) V ; 

The use of Mellin-Barnes integrals usually leads to rather large set of sums if compared to the 
sums needed to express the physics results. We refer to this method, since it uniquely allows to 
create the structures emerging. 6 All the Feynman-parameter integrals can now be performed 
and lead to a further factor of the kind 

r( ai Af + 6i(<t ) + fie) . . . T(a k N + b k (a a ) + r k e) 



T( Cl N + di(<7 ) + qie) . . . T(ciN + di{a a ) + q t e) 



(2.5) 

e z,fi,qj e Q 



6 A much more economic way consists in representing the different contributions to the Feynman parameter 
integrals in terms of generalized hypergeometric and related functions directly, cf. [6d] . 
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Here &«, di are linear functions of the Mellin-Barnes variables with integer coefficients. 

Now the k Mellin-Barnes integrals have to be performed. This can be done applying the 
residue theorem, after the according contours have to be properly fixed. In this way one obtains k 
or more infinite sums over the above rational function of T-functions, with according replacement 
of the sum-index in place of the variables a a . 

For these multiple sum expressions one may seek representations which are generalized hy- 
pergeometric series [7, 44] and generalizations thereof. Going to ever higher orders not all of 
these functions may be known by now but have to be introduced newly. They are naturally 
generated by the above integrals. 

At this stage one may wish to expand in e. If quantities without infrared singularities are 
considered, the maximal negative power in e is given by the loop-order of the diagram. Infrared 
singularities may enhance this order. The expansion has to be carried out up to the respective 
order required in the calculation. As only T-functions are expanded, which partly contain aiki, 
with ai,ki G N, single harmonic sums of the type Sb{aiki) emerge in products. The various 
infinite sums with index k\ nest these products, which finally leads to nested harmonic sums and 
multiple ^-values r their respective generalizations in higher orders, as described above. 7 The 
according sums can be uniquely found investigating the corresponding recurrences in A" over 
riE-fields, as possible with the code Sigma, [46]. 

The harmonic sums, which occur in the calculations of single scale quantities up to a given 
weight can be represented over a basis of functions, which is independent of the respective 
quantity being considered. In the following we identify these basis elements up to weight w = 5. 

3 Basic Definitions 

We will consider one-parameter real-valued functions x G ]0, 1[. Partly this set will be extended 
to special distributions f(x) G T>'(]0, 1[) [14] occurring in field theoretic calculations. Their 
Mellin transform is defined by Eq. (1.2), where AMs a positive, suitably large integer. According 
to the emergence of the respective Mellin transform in quantum field theoretic calculations, 
such as the local light cone expansion [39], A" is either an even or odd integer. Factors of (— 1)^, 
which emerge in the following, have therefore a definite meaning and are either equal to +1 or 
— 1. Both branches may be continued analytically in A" from N to C according to Carlson's 
theorem [47]. The distributions f(x) considered in the following are differentiable functions in 
the class C°°(]0, 1[) or ^-distributions 5^(1 — x), k > in T>'[0, 1] [14]. Furthermore, we 
consider functions with a +- prescription, 



with / G C°°(]0, 1[). The Mellin transform of the latter functions are harmonic sums in case 
the functions f(x) are harmonic polylogarithms [12] weighted by a factor l/(l±x), which we will 
consider throughout the present paper. The Mellin transform of the 5^(1 — ^-distributions 
yield 



7 First non-harmonic sums were found in individual diagrams in [45]. They canceled, however, in the final 
result in the MS-scheme within each color factor. 




(3.1) 




1=0 



(3.2) 
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In all cases below which are not of the form f(x) = ln m (x)/(x), or a linear combination of these 
terms, the function can be expanded into a Taylor series for x G [0, 1[. 



f( x ) = ^2akX , a k = — □ — , \/a k e R 



k=0 



k\ 



(3.3) 



The Mellin transform of f(x) is 

[ 1 dxx N f(x) = f:—^ 



+ k 



(3.4) 



All the above Mellin transforms are meromorphic functions in N. For the case f(x) = ln m (x)f(x) 



JO 



dx x N hi m (x)f(x) 



dx x f(x) 



(3.5) 



holds accordingly. 

The class of functions f(x) is formed out of polynomials of harmonic polylogarithms [12]. For 
physical applications in the massless case it will turn out that all contributing harmonic sums 
up to weight w = 5 can be represented by polynomials out of Mellin transforms of the type 



M 



x-l 



+ 



(AT), or M 



m 

x + l 



(N), 



(3.6) 



where f(x) is a polynomial of Nielsen integrals. In the analytic continuation one may consider 
Mellin transforms at rational arguments, like N/2 and study the mapping 



M 



f{x 2 ) 
1 — X 2 



(N) -> M 



1 — x 



N- 1 



(3.7) 



for /(l) = 0. The factor decomposition of the denominators yields 



1 - x 2 $i(x)$ 2 (^) 



1 1 
+ 



1 + x 



X 



(3.8) 



where ^\^{x) denote the first two cyclotomic polynomials [48] 8 . The representation of a har- 
monic polylogarithm f(x 2 ) in terms of other harmonic polylogarithms fi(x) may be worked out 
for non- negative indices, [12]. We will give the explicit relation for the case of the Nielsen 
integrals S n>p (x) [4] 



The following relation holds 



/ —W l - 1 (z)ln p (l-xz) . 



(n — l)\p\ Jo z 



(3.9) 



(n- l)\p\ ^ \IJ Jo z \xj K J 



(1 + z). 



(3.10) 



In massive calculations higher cyclotomic polynomials may occur as well. Examples are the fractional values 
of the /3-function (5.8), as /3((x+ 1)/2) = 2M[l/(l + z 2 )](i), /3((x + l)/3) = /3(x + 1) - M[(z - 2)/(z 2 - z + l)](x) 
for the 4th and 6th cyclotomic polynomial, cf. [27]. See also [49]. 
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For n — 1 one obtains 

-S liP (* 2 ) = 5 1|P (x) + 5 liP (-x) + E (J j / o — ln p -'(l - zx) ln'(l + **) . (3.11) 

Even in this simple case for Nielsen integrals, which are not classical polylogarithms, the mapping 
of the argument x — > x 2 does not give a closed relation in the same class of functions anymore. 
However, relations of this type may be used to express certain Mellin transforms, cf. [2], for 
examples up to w = 4. 



4 Classes of Harmonic Sums 

The number of all possible alternating and non-alternating harmonic sums of weight w, iV(w), 



is 

)W— 1 



7V(w) = 2 • 3 W ^ . (4.1) 

The number of basic sums in the case of all permutations can be calculated using the 1st Witt 
formula [17] for an alphabet of length I = 3 for w > 1, [15], 



7V basic (l) 

w rr \d 



7V basic (w) = - f^) 3 d , w > 2 . (4.2) 



' w 



The alternating and non-alternating harmonic sums [2,3] are generated by Mellin transforms of 
nested integrals based on the three w = 1 functions 

(4.3) 



1 + X x 1 — X 

Therefore, the alphabet out of which the index of the respective harmonic sums is formed as a 
word is of length I = 3. The harmonic sums form a sub-set w.r.t all possible index-sets over 
the alphabet { — 1, 0, 1}, since the letter cannot appear first. The counting relation is given by 
Eq. (4.1). 

Investigating the structure of the single-scale quantities in QCD and QED, cf. [6,7,10,29,31, 
50] to 2- and 3-loop order in more detail, it turns out that harmonic sums with index i = { — 1} 
do not contribute. One finds that the number of all harmonic sums of weight w, which do not 
contain any index % = —1, is obtained by expanding the following generating function [51] 

L "-'' E%d( w K. (4-4) 



1 -'' w ,1 



with 



= i [(1 - v/2)" + (1 + V2)™] = g (») 2' . (4.5) 

A r _,{_ 1 }(w) obeys the recursion relation 

N^ 1} (w) = 2 • A^ { _ 1} (w - 1) + A^ { _ 1} (w - 2) (4.6) 
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with the starting values N^{_iy(l) = 1, iV-,{_ 1 }(2) = 3 . Another interesting sub-class of the 
harmonic sums is that of the non-alternating harmonic sums. These sums are formed over an 
alphabet of length I = 2 [16]. Similarly to (4.2) one obtains 



Ni > (w) 

sic / 
v i > OV 



■>w— 1 



1 , 
1 

w 



5> 

d\w 



2 d , w > 2 . 



(4.7) 
(4.8) 



In Table 1 the number of sums is given in dependence of the weight w. Here the number 
of a-basic sums is the number of sums obtained using the algebraic relations [15]. 9 As shown 
in Table 1 the number of a-basic sums for all harmonic sums, Abasic(w), is larger or equal then 
the number of all harmonic sum, which do not contain { — 1} as an index. The corresponding 
number of basic sums is given by : 



i\^f- c i } (w) 



2 _ /w 



w 



a /-basic 
^{-1} 



(d). 



(4.9) 



/_ n(w) are partly connected to the number of a non-rooted tree studied in [52] for 



Values of N basic 

w = 2 ... 7, but obey a different counting rule otherwise. 

We mention that relation (4.2) applies also to other objects, as the harmonic polylogarithms, 
where the respective numbers up to w = 8 were determined empirically in [12]. Here we may 
apply it even for w = 1, unlike the case for harmonic sums. 



Weight 






Nu 


mber of 






Sums 


a-basic sums 


Sums ->{ — 1} 


a-basic sums 


Sums % > 


a-basic sums 


1 


2 


2 


1 


1 


1 


1 


2 


6 


3 


3 


2 


2 


1 


3 


18 


8 


7 


4 


4 


2 


4 


54 


18 


17 


7 


8 


3 


5 


162 


48 


41 


16 


16 


6 


6 


486 


116 


99 


30 


32 


9 


7 


1458 


312 


239 


68 


64 


18 


8 


4374 


810 


577 


140 


128 


30 


9 


13122 


2184 


1393 


308 


256 


56 


10 


39366 


5880 


3363 


664 


512 


99 



Table 1: Number of harmonic sums, number of sums, which do not contain the index { — 1}, number 
of sums with positive indices and the respective numbers of basic sums by which all sums can be 
expressed using the algebraic relations (a-basic sums) [15] in dependence of their weight. 

9 Later on we will derive a basis representation for the harmonic sums including the structural relations. Since 
in this representation we make a special choice, we list the corresponding algebraic relations w.r.t. to this basis 
for the dependent sums in Appendix A. 
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5 Sums of weight one 

The harmonic sums of weight w=l are generated by the Mellin transform of the first two cyclo- 
tomic polynomials [48] 10 (after regularization of the integral) 



1 -x 



and 



1 + x 



by 



M 



1 — x 



M 



1 + x 



(AO 
(AO 



-S 1 (N) = -^(N + l)- lE 
'-l) N [S- l (N) + ln(2)}=(3(N+l) . 



(5.1) 

(5.2) 
(5.3) 



Here ^ E denotes the Euler-Mascheroni constant. The representation of the integral in Eq. (5.2) 
in terms of the ^-function was first given by Legendre [53]. The Mellin transform for an 
integrable function x N f(x) obeys the identity 



M [/ (x a )} (AO 



M[f(x) 



N + 1 - a 



, a G R, a > . 



This relation modifies in case of f(x) being a +-distribution. For a = 2 one obtains 



M 



1 — x 



1 — X 2 
N -1 



1 
2 

= M 



1 1 

+ 



1 — x 1 + x 
1 



1-x 



(N) + M 



1 +x 



(AO + ln(2) 



(5.4) 

(5.5) 
(5.6) 



The latter relation can even be analytically continued in closed form since, see [2], 



-1> 



with 



(3(N) 



IE 



-^(N)- lE + (3(N) + \n(2) , 



N + 1 



cf. [27,54]. The duplication relation for the ^-function yields 



N+l 



+ ln(2) 



(5.7) 



(5.8) 



(5.9) 



Here fractional arguments in the ^-function emerge quite naturally. For k > 1 the harmonic 
sums S-k(N) and Sk(N) obey, respectively : 



S-k(N) = 
S k (N) = 



(fc-1)! 
(_l)(*-i) 

(fc-1)! 



(_l)^(*-i)(jV + l) - 



2 k- 



(5.10) 
(5.11) 
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These polynomials are the real-valued factors of x — 1 with N > l,JVsN. 
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with (k the Riemann ^-function. All these sums are expressed by the ^-function and its deriva- 
tives. We therefore conclude that due to the above rational-argument and differential relations 
there is only one basic single harmonic sum, Si(N). 

One may wonder whether differentiation relations can also be of importance for multiple zeta 
values. This is indeed the case [55]. The multiple zeta values can be generalized to multiple 
Hurwitz zeta values [56] defined by 



0= I Mfeft) • < 5 ' 12 > 



fci= 



iteratively with 



Here a £ C so that the above sums are defined. One may differentiate for a and perform 
then the limit a — > 0. Multiple Hurwitz zeta values obey stufhe relations, using the diction of 
Ref. [23]. Furthermore, one may form shuffle, duplication and more general algebraic relations 
among them, similar to the case studied for the multiple zeta values in Ref. [24]. Differentiating 
these relations builds classes of these relations which may be somewhat easier obtained as working 
out the relations directly for the multiple zeta values. It is not excluded that at very high weights 
even more relations between multiple zeta values can be found than obtained in [24], although 
for lower weights, i.e. w < 12 this is not expected. 

The differential quotient of a harmonic sum is always of the form 

-^Sol(lV) = S a - 2 (N) - Cg, lim S a - 2 (N) = ( a - 2 . (5.14) 

This is due to the fact that (5.14) can be represented by a Mellin transform which contains x N 
as a factor of the integrand being continuous in i £ ]0, 1]. The Mellin transform vanishes for 
TV -> oo. 

As an example we consider single harmonic sums. Differentiation yields 

jfiS k (N) = -k[S k+l (N)-C(k + l)} (5.15) 



_d_ 

In 



iys k (N) = -(-irk 



s_ (w) (iv) + (i-^)c(fc + i) 



(5.16) 



Differentiating the Hurwitz zeta value C a (^i) an d taking the limit a — * yields 



lim^C = -cC +(c+ i), OO (5.17) 
lim^C = cC-(c-i), c<0. (5.18) 



6 Sums of weight two 

Among the six harmonic sums four decompose into polynomials of single sums, which may be 
expressed through the tp^(z) functions as shown in the previous section. Two sums, S'_i ) i(iV) 
and Si t -i(N) remain. Due to the simplest shuffle relation by Euler [57] 

S a ,b(N) + S b , a (N) = S a (N)S b (N) + S aAb (N) (6.1) 
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with 



a A b = sign (a) sign (6) [\a\ + \b\] (6.2) 
these sums are connected, which implies the following relation for the Mellin transforms 



M 



ln(l — x) 



1 + x 



(TV) 



-M 



ln(l + x) 
1 +x 

N+l 



(N) 



+(-l) JV+i ^(^S^N) + S- 2 (N) 
+ [Si(N) - S_,(N)} ln(2) - ln 2 (2) + ^2 



(6.3) 



cf. [2]. Similar to (5.5) one may decompose 



ln(l — x 2 



1 f ln(l - x) ln(l - x) ln(l + x) ln(l + x) 



X 



which yields 



M 



ln(l — x) 
1 — x 



N-l 



= M 



+ M 



1 + x 



dn(l — x) 
1-x 

ln(l + xj 
1 + x 



X 



1 + x 



(6.4) 



(iV) + M 
(N) + M 



ln(l — x) 
1 +x 

\n(l + x)\ 
1-x ) 



(N) 
(AO 



-2C2 + ln 2 (2). 



(6.5) 



Since the l.h.s. of (6.5) and the first and fourth term in the r.h.s. are polynomials of single sums, 
the remainder two terms are related. However, (6.4) does not yield a new relation. Therefore 



Fi(N) := M 



\n(l + x) 
1 + x 



(AO 



(6.6) 



is the first non-trivial Mellin transform beyond the single harmonic sum Si(N). Although Fi(N) 
expresses a harmonic sum containing { — 1} as index, which belongs to a class being absent in 
the final results in physics applications, it may be of use to derive compact expressions in other 
cases later. 

As a historic remark, we mention that Nielsen [27] considered these functions as well using 
the notation f(jV), rj(N), ^(N) and f 2 (iV), 



rj(N) 
UN) 

-6 (AO 



M 



'ln(l — x) 



x-1 ) 
S ltl (N-l) 

~\n(l + x) - ln(2) 



(N-l) 



^'(N)-C2-^{N)+ lE f 



M 
M 
M 



x - 1 
ln(l + x) 

x + 1 
ln(l — x) 

x + 1 



(N-l) 



(N-l) 
(N-l) 



(6.7) 
(6.8) 

(6.9) 

(6.10) 
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They obey the relations 



z) + lE ] [4>(l-z)+ lE ] 


= 2C2-Z(z)-Z(l-z) 


(6.11) 


P(z) [i>(l-z)+ lE ] 


= -/5(l-z)ln(2)-6W-6(l-^) 


(6.12) 


P(z)[^(z)+ lE ] 


= ^) + ^)ln(2) + Ieg)-^(^ 1 ) 


(6.13) 


P(z)+P(l-z) 


7T 


(6.14) 


sin(7T2;) 


(3(z)(3(l - z) 


= rj(z) + 77(1 — z) 


(6.15) 




= rl>'(z) - 2r)(z) , 


(6.16) 



for z E ]0, 1[. Most of them were given in Ref. [27]. They result from algebraic relations between 
harmonic sums and corresponding integral representations. Let us illustrate this for (6.12). One 
may represent 

P(N + 1) [i>(N + 1) + lE ] = (-1) N S 1 (N) [S-!(N) + ln(2)] , (6.17) 

cf. [2]. Using Euler's relation (6.1) one obtains 

P(N + 1) ty(N + 1) + lE ] = (-If [S^N) + S- ltl (N) - S- 2 (N) + H2)S 1 (N)] , (6.18) 

which can be expressed referring to the functions P(N), P'(N), £i(iV) and ^(N). Relation (6.13) 
is a consequence of the properties of the logarithm. 11 



7 Sums of weight three 



There are eighteen harmonic sums of weight w = 3, [2]. Twelve of these sums are related to the 
remaining six harmonic sums by shuffle relations [15]. The latter are represented by the Mellin 
transforms 



M 



'\n 2 (l + xY 
1 -x 



(iV) 



M 



r ln 2 (l + x) 
1 + x 



(TV) 



M 



Li 2 (-x) 
1 — x 



+■ 



(AO 



M 



'LisOr) 
1 — x 



(AO 



(7.1) 



M 



Li 2 (-a:) 
1 + x 



(AO 



M 



Li 2 (x) 
1 + x 



(AO 



The first two functions in (7.1) correspond to the harmonic sums of the type S'_i i i ) _i(A r ) and 
Si, 1,-1 (A"), which are not connected algebraically. Using the general relation 



M[\n l (x)f(x)] (JV) = [/(*)] (N), 



(7.2) 



11 We correct misprints in Ref. [27], §75, third equation. The minus sign in front of j3'(x) has to be removed. 
The equation below Eq. (8) has to be replaced by (6.12). 
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the Mellin transform M [ln(x) ln(l +x)/(l + x)] (TV) can be calculated from Fi(N). We exploit 
Euler's relation for the sums Si - 2 (N) and S_ 2jl (N), [2], to express M [Li 2 (-x)/(l + x)] (iV) in 
terms of M [Li 2 (x)/{1 + x)] (N), 



M 



Li 2 (-x) 
1 + x 



(TV) = (_i)^|ln(2)[5 2 (7V)-5_ 2 (7V)] + ^-iW + 7C3 
+S. 1 (N)S 2 (N) + S_ 3 (TV) 1 + M 



2 ~~-xv--, ■ 4 
Li 2 (x) + ln(x) ln(l — x) 



1 + x 



(TV) 



Furthermore, the identity 

1 LiJx 2 ) 



LiJx) LiJx) LiJ-x) LiJ-x) 



2 n-2 1 _ x i 

holds and one obtains 

'Li 2 (-x)\ 



1 — x 1 + x 



M 



x- 1 



(TV) 



1m 

2 

-M 



'Li 2 (x) 
x- 1 

' Li 2 (x) 
, x — 1 



1 — x 



TV- 1 



(JV) + M 



1 + x 



(7.3) 



(7.4) 



+ M 

Li 2 (x 



x + 1 



Li 2 (-x) 
x + 1 

(TV) 



(TV) 



i C 3 + ^C 2 ln(2). 



(7.5) 



Eq. (7.5) can be used to express M (Li 2 (— x)/(l — x)) + (TV) in terms of the remaining Mellin 
transforms. 

Four basic functions emerge at weight w = 3 



F 2 (N) 
F 3 (N) 
F*(N) 
F 5 (N) 



M 
M 
M 
M 



ln 2 (l + x) -ln 2 (2) 



1 -x 
r ln 2 (l +x 



(AO 



1 + x 
Li 2 (x)" 
1 — x 



+ 



(AO 
(TV) 



Li 2 (x) 
1 +x 



(TV). 



(7.6) 
(7.7) 
(7.8) 
(7.9) 



Here we have chosen a form for F 2 (N) which yields a particular simple asymptotic representation 
as A" G C, | TV | — > oo. The basic sums at weight w = 3 which do not contain the index { — 1} 
are thus given by 



5_ 2 ,! (TV) 



1) N+1 F 5 (N) + C 2 S_!(TV) - ^Cs + C 2 m(2) 



5 



8 



S 2 ,i(N) = F 4 (TV)+C 2 5i(iV) . 



(7.10) 
(7.11) 
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8 Sums of weight four 



Out of the 54 harmonic sums which emerge at weight w = 4, 38 harmonic sums can be ex- 
pressed by 16 sums through shuffle relations [15]. The latter sums are related to one Mellin 
transform for each index set ( — 1, 1, 1, 1), ( — 1, — 1, — 1, 1), ( — 1, — 1, 1, 1), ( — 2, 2), ( — 3, 1), (3, 1), 
(-3, -1), (3, -1), (-1, -1, -2), (-1, -1, 2), (1, 1, -2) and (1, 1, 2) and to two Mellin transforms 
for the sets (—1, 1, —2) and (—1, 1,2). Various of these harmonic sums contain an index { — 1}. 
Since sums of this type do not occur in the known physics applications we will discuss sums 
without this index only from here on. After the algebraic relations are exploited, the seven sums 
with index sets (2, 1, 1), (-2, 1, 1), (3, 1), (-3, 1), (-2, 2), (4), and (-4) remain at this weight. 

Unlike the case for weight w = 2, the harmonic sums with symmetric index pattern at weight 
w = 4, S-2,-2(N) and S^A) relate two of the above Mellin transforms. They have the 
representation, [2]. 



S-2,-2( N ) = 



S 2 , 2 (A) 



'\n(x)U 2 (-x) -2U 3 (-x) 
x-1 



M 

+|[5 2 (JV)-5_ 2 (JV)]-|cs5i(JV) 



(A) 



M 



P'(N + 1) 



\N ■ 



<2 



/ ln(x)Li 2 (a;) — 2Li 3 (a;) 
x-1 



2 -^ (3) (iv+i) + jkf 



(N) + 2( 3 S 1 (N) 



= M^'(iv + i)-C2] 2 -^ (3) (iv + i) + ?C 2 2 



6 



(8.2) 



Since the Mellin transform of the distribution [Li 2 (— x)/(x — 1)]+, (7.5), depends on Fi t2 ^(N) 
and M (Li 2 (x)/(x — 1)), (A) is a basic function their derivative w.r.t. A is known and the 
Mellin transforms of [Li 3 (—x)/(x — 1)] + and [Li 3 (x)/(x — 1)] + can be expressed by (8.1,8.2). 
One may use the relation 



S 3 ,_i(A) + S- 1>3 (N) = S^(N)S 3 (N) + S_ 4 (A) 



(8.3) 



to express M [Li 3 (— x)/(l 
Li 3 (-x) 



M 



1 +x 



(AO = 



x)} (N) 




(A) + M 



ln(x)U 2 (x) + ln 2 (x) ln(l - x)/2 
1 +x 



+(-!)" <nQ - oCsln(2) + S^(N)S 3 (N) + S. 4 (N) 



ln(2) [S_ 3 (N) - S 3 (N)} + ^C 2 S- 2 (A) - ^Cs^-i(A) 



(A) 



-(5_ 1 (A)+ln(2))C 3 • 



(8.4) 
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Furthermore, the relation 



M 



Li 3 (a;) 
1 — x 



M 

+ M 

29 2 
+ 40 C2 



'U 3 (x) 



1 — x J 

\ l-x 
2Li 4 (- 



(N) + M 



U 3 (x) 



(N) + M 



1 + x 
Li3 



(N) 



(8.5) 



1 



x 



(TV) 



^ln(2)C3 + ic 2 ln 2 (2)-lln 4 (2) 



holds. However, similar to (6.5), it does not lead to a further reduction since 
M [Li 3 (-x)/(l + x)] (N) and M [Li 3 (x)/(1 + x)] (N) enter in the same way as in (8.4). 

Square argument relations may be considered for more general Nielsen integrals as well. For 
Si, 2(^)/(l ~~ x ) one obtains the relation 



M 



'Sio(x) 



X 



N - 1 



= M 
+ M 
+ M 



+ 



Jo 



' Sl,2{ X ) \ 

l-x) + _ 

Si, 2 (-^) 
1 — X 

d x S ^ 



(N) + M 

(iV) + M 
(N) + M 



Si 2 (x) 



1 + x j 
Si 2 (-x 



(N) 



1 + x 



(N) 



1 + x, 



(AT) 



1 + x 

which refers to the functions Si i 2(±x)/(1 ± x) and h(x)/ (1 ± x). h(x) is given by, cf. [30], 

1 



(8.6) 



h(x) 



-SiAx ) - SY 2 (x) - Si 



-x 



(8.7) 



The function 7 1 (x)/(l — x) does not occur up to w = 6 for single scale quantities. Therefore, we 
do not consider square argument relations beyond those for Li n (x) in the following. 
At weight w = 4 the following basic functions contribute : 



F 6b (N) 
F 7 (N) 
F 8 (N) 



:= M 

:= M 

:= M 

:= M 



Li 3 Qr) 
1 + x 

( SiM 

\x-l 

Sl,2(x) 

1 + x 



(N) 



(AT) 



(TV) 



(AT). 



(8.9) 
(8.10) 
(8.11) 



The number of basic Mellin transform occurring in physical processes turns out to be even 
smaller. In the case of the one-loop anomalous dimensions and the Wilson coefficients in first 
order in the coupling constant, at most single harmonic sums contribute. In the case of the two- 
loop anomalous dimensions F^(N) emerges as the first non-trivial function, cf. [1]. The case 
of the 2-loop coefficient functions has been analyzed systematically in Refs. [10] for all Wilson 
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coefficients in polarized and unpolarized deeply inelastic scattering, the polarized and unpolar- 
ized Drell-Yan process and the fragmentation functions, as well as the coefficient functions for 
hadronic Higgs- and pseudo-scalar Higgs-pro duct ion. In particular Mellin transforms being 
associated to harmonic sums with an index set formed of chains of the elements — 1 and +1 do 
not contribute at this level. This does not exclude that functions of this type do not emerge in 
higher orders. 

9 Sums of weight five 

We will consider the sums of different depth separately and assume that the algebraic relations 
have been exploited already. 



9.1 Twofold Sums 

The following sums contribute : S±4,i(iV), S± 3: ± 2 (N). We obtain the representations 



S 4A (N) = -M 



x — 1 



(AO + S^NX* - S 2 (N)( 3 + S 3 (N)( 2 



S_4,i(A0 = (-1)^ +1 M 



Li^x) 
x + 1 



(AO + S-!(A0C4 - S- 2 (N)( 3 + S_ 3 (A0C 2 



+C 4 ln(2) + -C 2 C 3 -^C5 



S 2 ,s( N ) 



M 



'ln(x) [5 1)2 (1 - x) - Ca] + 3 [5 1)3 (1 - x) - Q 



x - 1 



(N) 



(9.1) 



(9.2) 



(9.3) 



S-2, 3 ( N ) = (— 1) M 



\n(x) [S h2 (l -x)- Cs\ + 3 [S 1[3 (l -x)- Q 
1 +x 



(AT) 



+3&S-i(JV) + |c 5 + 3C 4 ln(2) - ^Cs 



S2,-s(N) = (-1) 



N+l 



M 



l + x 



■ ln 2 (:r)Li 2 (-x) - 2 ln(x)U 3 (-x) + 3Li 4 (-x) 



(9.4) 



(9.5) 



(AT) 



+?C3 [S. 2 (N) - S 2 (N)\ - yC4^-i(iV) - ^C 5 - yC 4 ln(2) + C2C3 (9.6) 



S- 2 ,s( N ) 



-M 



1 (\ 



x-\\2 



ln 2 (x)Li 2 (— x) — 2 \n(x)Li 3 (—x) + 3Li 4 (— x] 



(N) 



+ JjC3 [S 2 (N) - S. 2 (N)] - yC 4 Si(A0 



In the above integrands we use the relations 

Si l2 (l-a;) = -U 3 (x) + ln(x)Li 2 (a;) + ^ln(l -x)\n\x) + ( 3 



(9.7) 



(9.8) 



51,3(1 -x) = -Li 4 (x) + \n{x)U 3 {x) - ± ln 2 (:r)Li 2 (:r) - - ln 3 (x) ln(l - x) + (a , (9.9) 
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which reduce the representation. Since all Mellin transforms 



Lifc(=hr) 
x± 1 



, Jfe < 3, 



are known, the only new ones are those of the functions 

Li 4 (±x) 
x±l 

and the other terms follow by differentiation. Due to 



(9.10) 



(9.11) 



S^N) = (-1) N M 



Li 4 (-ar) - ln(x)U 3 (-x) + ln 2 (a;)Li2(-a;)/2 + ln 3 (x) ln(l + x)/6 

x+ 1 



7 1 7 29 

+ -C4 [S-i(N) - S 1 (N)] - -C 2 C 3 + gC 4 ln(2) + -Cb 



M 



Si,s(l-x) 



x-1 

Eqs. (9.2, 9.12) allow to express 



+ §Cl^i W + C2C3 - 2Cs 



(AO 
(9.12) 
(9.13) 



M 



Li 4 (— x) 
1 +x 



(N) 



(9.14) 



but (9.13) does not lead to a new relation. 

In some of the harmonic sums Mellin transforms of the functions 

Lifc(-rr) 
x±l 



(9.15) 



emerge. For odd values of k = 2/ + 1 the harmonic sums S 1 -^-i)(N), S^k-i),i(N) and S-i-i(N) 
allow to substitute the Mellin transforms of these functions in terms of Mellin transforms of basic 
functions and derivatives thereof. For even values of k this argument applies to M[Lifc(— x) /(l + 
x)](N) but not to M[Lifc(— x)/(l +x)](N). In the latter case one may use relation (7.4). Since 
in massless quantum field-theoretic calculations both denominators occur, one may use this 
decomposition based on the first two cyclotomic polynomials, cf. [48], and the decomposition 
relation for Li fc (a; 2 ). The corresponding Mellin transforms also require half-integer arguments. 
The relation 



2*-i 



M 



'Li fc (x 2 ) 
x 1 - 1 



N- 1 



M 

+ M 



L 



Lifc(rr) \ 
x-1 J 

Lifc(-ff) ' 
x- 1 

dx^-l 

1 + X 



(N) + M 



'Li fe (a;) 



+■ 



(N) + M 



x + 1 / + 
Life (-a;) 
x + 1 



(N) 



(N) 
(9.16) 



determines M[Life(— x)/(l + x)](N). For k — 4 the last integral in (9.16) is given by 



Jo 1 + x 



2 2 s ) 
- ln(2)C 2 2 + 3C2C3 - T C 5 



(9.17) 
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and one obtains 

Li 4 (— x) 



M 



x+1 



(TV) = 



1 



8 

+ M 
1 



M 



( Li A (x) 
\x-l 
Li 4 (— x) 
x — 1 



N — 1 



(N) - M 



+ M 

Li 4 (a;) 
x + 1 



Li 4 (x) 
x — 1 

(AO 



(N) 



^ln(2)-^C3 + |c 5 . 



The new functions are : 



F 9 (N) 
F 10 (N) 



M 



= M 



/ Li 4 (rr) 
\x-l 
Li 4 (x) 
1 + x 



(N) 



(N). 



(9.18) 

(9.19) 
(9.20) 



9.2 Threefold Sums 

At this depth the sums £±3,1, i( AO, £2,2,1(^0, £-2,-2,1 (^0, £-2,2,1 (^0 and S 2 -2,i(N) contribute. 
One derives the following representations. 



S-3,i,i(A0 



M 



'£2,2(2;) 



(-l^M 



x — 1 J 

S 2 o(x 



l + x 



(N) + ^S 2 {N)-^S l {N) 
(N) + C 3 S-2(N)-^S^(N) 



(9.21) 



£2,2,1 (N) 
£-2,1,-2^) 



+|CsC2 - ^C 4 ln(2) - ^Csln 2 (2) + ^ 2 hi 3 (2) 
+ | C ,- 21n(2) Li 4 Q)-2Li 5 (i)-lln 5( 2) 

(N)+(2S2AN)-^-ClS 1 (N) 



-M 



M 

1 

~2 



^(^-Li 2 .^)^ 

A 1 

Hx)S lt2 (-x)-Ul(-x)/2 ' 
x - 1 



(9.22) 
(9.23) 



(AO 



C2 [£-2,i(iV) - S-2,-i(A0] - [\Cs - \ ln(2)C 2 ] [S. 2 (N) - S 2 (N)\ 



S. 



+gC 2 2 £i(iV) 

(—1) N+1 M 



(9.24) 



2S 2 , 2 (x)-L%(x)/2 
1 +x 



(iv) + c 2 £- 2 ,i(iv)--c 2 2 £-i(iv) 



+4Li 5 (- J + 4Li 5 (- j ln(2) - -( 5 - -C2C3 + ^ ln 5 (2) - -( 2 ln 3 (2) 



+i C 3 ln 2 (2)--C 2 2 ln(2) 



(9.25) 
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S 2 ,i,- 2 (AO 



1)*M 



ln(x)S lt2 (-x)-Uj(-x)/2 
1 +x 



(TV) 



-\C2 [S 2 ,m - S 2 ^(N)] - [|Cs - ^C 2 ln(2) 

1 1 77 1 1 1 

+ gC^-i(iV) - -^C 5 + yC 2 C 3 + gC 2 2 ln(2) . 



[S 2 (N) - S_ 2 (N)} 



(9.26) 



Here also the sums S± 2t _i(N) emerge, which is a consequence of the representation chosen. The 
associated Mellin transforms are accounted for, however, by the w = 3 basic functions, cf. [2]. 
In the threefold sums the basic functions 



emerge. 



Fn(N) 
F 12 (N) 
F 13 (N) 
Fu(N) 
F 15 (N) 
F 16 (N) 



= M 



M 



M 



M 



M 



M 



S2M 

\x-l 

S 2)2 {x) 
1 + x 
/Li|(x) 

U-i 



(TV) 



(AO 



(AO 



Li 2 Jx) 



(AT) 



Hx)S lt2 (-x)-Ul(-x)/2 + Q/8 



x-l 

2< 



(AT) 



\n(x)S lj2 (-x) -Li 2 2 (-x)/2 
1 + x 



(N) 



9.3 Fourfold Sums 

Here the two sums £±2,1,1,1 (AO contribute. They obey the representation 



S 2 ,i,i,i(N) = -M 



x-l 



+ 



(9.27) 
(9.28) 
(9.29) 
(9.30) 
(9.31) 
(9.32) 



(9.33) 



S- 2 ,i,i,iW 



1)^ +1 M 



Sl,3( X ) 

1 + x 



(N) + US^(N) + C 4 ln(2) - ^GCs - \t2 ln 3 (2) 



1 ^ - + ^G) + h in5(2) + Li5 G 



The final basic functions contributing up to weight w = 5 are thus 

'Si, 3 (xy 



F 17 (N) = M 
F 18 (N) = M 



\x-l 
1 + x 



(AT) 



(AT) 



(9.34) 

(9.35) 
(9.36) 
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Let us summarize the basic functions, which contribute up to the level of the 3-loop anoma- 
lous dimensions, which are w = 5 quantities. They are given by 



w = 1 : l/(x- 1) 

w = 2 : ln(l + x)/(x + 1) 



(9.37) 
(9.38) 
(9.39) 
(9.40) 



w = 3 : U 2 (x)/{x± 1) 
w = 4: Li 3 (x)/(x + l), 
w = 5 : Li 4 (x) /{x ± 1), 



5 li2 (x)/(x±l) 

S lj3 (x)/(x + 1), ^(^/(xil), 
[ln(x)^, 2 (-x)-Li2(-x)/2]/(x±l) . 



Lil(x)/(a; + l), 



(9.41) 



The five-fold sum S^i^i^A") decomposes algebraically into a polynomial out of single sums, 
[2]. The representation of the regularized Mellin transforms of the basic functions for complex 
argument is given in Appendix B. 



We investigated the relations between harmonic sums, which are related to their functional value. 
They extend the set of algebraic relations being implied by their quasi-shuffie property related 
to the index pattern. In quantum field theoretic calculations of massless single scale quantities 
up to 3-loop order and for massive quantities to 2-loop order only harmonic sums with indices 
i 7^ — 1 occur. We derived the structural relations for this sub- algebra up to weight w = 5. 
These are fractional argument-, integration-by-parts- and differentiation relations. The set of 69 
harmonic sums is reduced to 30 sums by the algebraic relations. The structural relations imply 
a further reduction to 15 basic harmonic sums. In physical applications these sums have to be 
known for complex arguments. We showed that these functions can be analytically continued to 
meromorphic functions with poles at the non-positive integers referring to their representation 
in terms of Mellin-integrals. The basic sums obey recurrence relations for complex arguments 
given by corresponding harmonic sums of lower weight. Separating divergent contributions for 
\N\ — * oo, oc Si(N), k > 0, k G N algebraically, the asymptotic representation of the harmonic 
sums was derived in analytic form. The recursion relations and asymptotic representations allow 
to calculate the basic functions for complex arguments analytically. In the physical applications 
to w = 5, 15 basic functions contribute, as the example of the 3-loop anomalous dimensions 
shows. This set is extended to 37 functions at w = 6 in case of the 3-loop Wilson coefficients, 
cf. [37]. Here, 20 functions emerge for w = 6 and two w = 5 functions M[(Li2(x)/(x — l))+](iV) 
and Nl[(Si t3 (x) / (x — 1)) + ](N) do also occur. They do not contribute to the 3-loop anomalous 
dimensions. We included their representation in the present paper. The investigation of a large 
number of massless and massive 2-loop problems showed that the basic functions are universal 
up to isomorphics. This is expected also for the 3-loop case. 
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A Algebraic Relations 



In this Appendix we list the algebraic relations w.r.t. the basis derived in the present paper. 
As the two-fold sums obey the simple relation (6.1) we list the sums of higher depth only. In 
the following we give the explicit representation of the harmonic sums without an index { — 1}, 
which contribute in physical single-scale processes up to w = 5. The single sums are given in 
(5.10, 5.11). In the following we drop the common argument N. 

Weight 2 Sums : 



Si, 



sl + s 2 



(A.l) 



Weight 3 Sums 



Si. 



(A.2) 



Weight 4 Sums 



5*1,-2,1 — — 25'_2,l,l + SiS-2,1 + S-3 : \ + S_2,2 

5*1,1,-2 — S -2,1,1 — SiS-2,1 — 5-3,1 — 5_2,2 + 5i5_3 + 5_4 + 2^- 2 ^1 + ^2 



51.2.1 — ~2S , 2,1,1 + 5i52,l + 53 ; i + 52,2 

51.1.2 — 52,1,1 — SiS 2t l — 53,1 — 52,2 + S3S1 + 54 + -S2 Si + 52 

1 1 2 1 1 2 1 4 

5i,i,i,i = -^S 4 + -S 2 + -S 1 S 3 + -S 1 S2 + —S 1 



(A.3) 

(A.4) 
(A.5) 

(A.6) 
(A.7) 



Weight 5 Sums 












"S'l-3,1 = 


— 25_3,i,i + 5i5_3,i + 5_4,i + 5_3,2 


+ ^5-3 








5l,l,-3 = 


5-3,1,1 — 5i5_3,i — 5_4,i — 5_3,2 + 5i5_4 + 5_5 




-s 2 




5l,3,l = 


— 25*3,1, 1 + 5i53,i + 54,1 + 53,2 










5l,l,3 = 


53,1,1 — 5i53,i — 54,1 — 53,2 + 5i54 + 55 + -S3 


si + S 2 








5-2,-2,1 = 


- [ — 5_2,l,-2 + 5_25_2,1 + 54,i + 5_2,-3] 










"S'l-2-2 = 


— [ — 5-2,1,-2 + S-2S1-2 + 5_3,-2 + 5i,4] 










52,-2,1 = 


— 52,1,-2 — 5_2,2,1 + 53,_2 — 5i,_4 — 5_25i,2 + 5^2,-2 + 52,-3 


+ 5i5_ 


5-2,1,2 = 


525-2,1 + 5_4 I + 52,1,-2 — 53,-2 + Si 4 + S-2S12 — S1S2 _2 


_ 5 2 ,_3 



(A.8) 
(A.9) 
(A.10) 
(A.ll) 

(A.12) 

(A.13) 

5-2,3 

(A.14) 

'-2,2 
(A.15) 
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Si, 



-2,2 



— 525-2,1 — £-4,1 — 5*2,1,-2 — £-2,2,1 + 5*3,-2 

—Sl,-4 — £-2£l,2 + S\S2,-2 + 52,-3 + 25i5_2,2 + 5_2,3 + 5_3,2 



5l,2,-2 — 5_25i ; 2 + £l,-4 — 5i5. 



2,2 



-2,3 



+ 5_ 



2,2,1 



— 25*2,2, 1 + 525*2,1 + 54,1 + 52,3 



52,1,2 
5l,2,2 

5l, -2,1,1 — — 35-2,1,1,1 + 5i5_2,l,l + 5_3,i,i + 5_2,2,1 + 5-2,1,2 

5l, 1,-2,1 — 35-2,1,1,1 — 5i5_2,l,l — 5_3,i,i — 5_2,2,1 — 5_2,1,2 



52,2,1 + — [525i,2 + 53,2 + 5i,4 — 5252,1 — 54,1 — 52,3] 



+ 2 [5l5i,_2,l + 52,-2,1 + 5l,-3,l + 5i,_2,2] 



5n 



— 5_2,i,i,i + -5i I 5i,i,_2 + 5_2,i,i — -£1,-2,1 



+ 2 (52,1,-2 + 5l,2,-2 + 5l,l,-3 + 5_3,i,i + 5_2,2,1 + 5_2,1,2) 



6 



(52,-2,1 + 5l,-3,l + 5i,_2,2) 



5l,2,l,l — — 352,1,1,1 + 5i52,l,l + 53,1,1 + 52,2,1 + 52,1,2 



5i, 1,2,1 — 352,i,i,i + 5i 



2~5i,2,i - 5 2 ,i,i 



(A.16) 
(A.17) 
(A.18) 

(A.19) 
(A.20) 

(A.21) 



(A.22) 
(A.23) 



+ 2 [£1,3,1 + £1,2,2 - £2,2,1] - £3,1,1 - £2,1,2 



£11 



1 



1,2 



1 



-£2,1,1,1 + g£l [£l,l,2 + £2,1,1 - 2^1,2,1 J + 3 

+ -£l,2,2 + 7^£2,2,1 + £l,l,3 — 7T£l,3,l + £3,1,1 



252,i, 



(A.24) 

(A.25) 
(A.26) 



Some of the sums of lower complexity can be further expressed using algebraic relations. 
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B The Basic Functions for Complex Arguments 



To derive the analytic representation of the Mellin transforms of the basic functions we first 
express them in terms of factorial series and a rest part which can be decomposed algebraically. 12 . 
We show that both parts are meromorphic functions with poles at the non-positive integers. 
They obey recursion relations in which functions of the same class but lower complexity occur. 
Finally we derive the associated asymptotic representations in analytic form. To make use of 
this formalism the basic functions have to be mapped. The corresponding functions lay in the 
same equivalence class. First we remind some properties of factorial series. Then we define 
the functions associated to the basic Mellin transforms and derive their recursion relations and 
asymptotic representations. 



B.l Asymptotic Representations of Factorial Series 

The Mellin transforms 

Q(z) = f 1 dt ip(t) t z - 1 (B.l) 
Jo 

with tp{t) one of the above basic functions, suitably transformed, is an integral representation of 
a factorial series Q(z) [27,28] in the complex plane z G C. The basic functions are mapped 
such, that (p(t) can be expanded into a Taylor series around t — 1, 

oo 

p(l -*) = !>*** ■ (B.2) 

k=0 

As will be shown below, the remainder terms can be expressed in terms of known harmonic sums 
of lower complexity. 

For Re(z) > 0, Q(z) is given by the factorial series 

" w = + d ■■'■(« + *) • (a3) 

fl(z) has poles at the non-positive integers and one may continue Q(z) analytically to values of 
z G C as a meromorphic function. Let us now derive the asymptotic representation 

oo 

m - e - k ■ (b.4) 

k=l z 

As shown in [27] the coefficients are obtained by 



di = ai 

k—2 

a k = ]T(-l)^a fc _,, (B.5) 

1=0 



where C l k are the Stirling numbers of 2nd kind, 

1 k-2 



12 In [58] asymptotic relations for non-alternating harmonic sums to low orders in 1/N k were derived. Our 
algorithm given below is free of these restrictions. The main ideas were presented in January 2004 [59], see 
also [29]. 
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e k = i (b.7) 

4 = € l k+1 -k€[- + \. (B.8) 

Not for all harmonic sums the asymptotic representation can be found in the way outlined 
above using the associated Mellin transform [2, 15]. In the case of Si(N) the function ip(N + 1) 
contains a logarithmic contribution 

00 z 



1 ^ B 



Here denote the Bernoulli numbers, which can be obtained from the generating function [60] 

T oo k 

Also ip(z) and its derivatives are meromorphic functions with poles at the non-positive integers. 
One may show that nested harmonic sums with an index pattern 

{1,1,... ,l,a,...}, \a\>l (B.ll) 

contain a term oc S±(N) as \N\ oo, where k is the number of l's left to a. These terms 
can always be separated algebraically and result in a groth oc \n h (z) similar to the behaviour in 
(B.9). After this separation was performed the remainder terms can be expressed by factorial 
series. 

The function v{z) is associated to ip(z) by 

u(z) = ln(z) - rj)(z) (B.12) 
and removes the logarithm in (B.9). It obeys the following factorial series [27] 

where ip n (x) denotes the Stirling polynomial which obeys the recursion 

(a; + 2)ijj n (x + 1) = (x - n)^ n (x) + (x + l)^n-i(^) , (B.14) 
with the initial values 

Mx) = \, ^i(x) = i(3a; + 2) . (B.15) 

The function (3(N) associated to the harmonic sum S-i(N — 1) is represented by the factorial 
series 

oo K| I oo (_-\\k 
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To derive asymptotic representations for the basic function Fi(z)\i=i..is some of the functions 
have to be mapped into a form, which is regular at z — 1, which is indicated by a hat below. 
We consider the following Mellin transforms : 



F 2 
F 3 
F 4 

h 

Fe a 
Feb 

A 

F 8 
F 9 
Fio 
F n 

F\2 
F13 

F u 
F 15 

Fn 



M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 



ln(l + x) 



1 + x 
ln 2 (l + :r) -ln 2 (2) 

l-x 
ln 2 (l + 1 



1 + x 
Li 2 (l-rr) 

1 — x 
Li 2 (l-x) 

1 + x 

1 — X 

1 + x 
U 3 (l-x) 

1 — X 
Li 3 (l-a;) 



1 + X 




SiflO- - 


x) 


l-x 




SiflO- - 


x) 


1 + x 




s 2 , 2 (i - 


x) 


l-x 






x) 



[z) 



1 +x 
Li^(l-rr) 
1 — x 



(z) 
(z) 
(z) 
(z) 
(z) 



Ul(l-x) 
1 + x 

Hx)S 1 , 2 (-x)-Uj(-x)/2 + Cj/8 
1 — X 

\n(x)S h2 (-x) - Li 2 2 (-x)/2~ 



(z) 



Li 4 (l — x) 
l-x 



1 + x 



(z) 



(z) 
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Fis(z) 



M 



Li 4 (l — x) 
l + x 



(B.35) 



B.2 Representation of the Basic Functions 

The numerator functions are related to the original basic functions, cf. [61], by 
Li 2 (x) = — Li 2 (l — x) — ln(x) ln(l — x) + ( 2 



Li 3 (x) 
Li 4 (x) 



-51,2(1 - x) - In (x) Li 2 (l -x)-- ln 2 (x) In (1 - x) + ( 2 In (x) + ( 3 



-Li 3 (l -x) + ln(l - x)Li 2 (l - x) + - ln 2 (l - x) \n(x) + Cs 
-S 1>3 (1 -x)- Hx)S lt2 (l -x)- 1 - ln 2 (x)Li 2 (l - x) 
-- \v?(x) ln(l -x) + ^ ln 2 (x)C 2 + ln(x)C 3 + C 4 



(B.36) 

(B.37) 
(B.38) 

(B.39) 



U 2 2 (x) = 
Sv*{x) = 



-Li 4 (l -x) + ln(l - x)Li 3 (l -x)-- ln 2 (l - x)Li 2 (l - x) 

-iln 3 (l -x) ln(x) + (4 
6 

Li 2 (l -x) + 2Li 2 (l - x) ln(l - x) ln(x) - 2Li 2 (l - x)( 2 
+ ln 2 (l - x) \n 2 (x) - 2 ln(l - x) \n(x)( 2 + C2 

-5 2 , 2 (1 -x) + ln(l - x)5 li2 (l - x) + I 



(B.40) 
(B.41) 



- [Li 3 (l -x)- ln(l - x)Li 2 (l -x)- Cs] ln(x) + - ln 2 (l - x) ln 2 (x) . (B.42) 

In the above we arranged the representations such, that the remainder Mellin transforms can be 
obtained either by algebraic relations or differentiation of Mellin transforms being derived above 
or are known otherwise, cf. [2,31,62]. We remind the following relations : 



M[ln(l - x)Li 2 (l -x)](N) 



M[\n 2 (l - x)U 2 (l - x)}(N) 



M[ln(l -ar)Li 3 (l - x)](N) = 



M[ki(l-x)S 1>2 (l-x)](N) = 



N + l 

-ls 2tl (N + l) 



S^N + l)S 2 (N + 1) + S 3 (N + 1) - C 2 5i(iV + 1) 

(B.43) 

6.$'ju. 2 (A' L) r I- , V-l)-25 2! u(.V) 

(B.44) 



N + l 

-6( 2 S 1>1 (N + 1)-2( 3 S 1 (N + 1) 
1 



N + l 



35 1)1 , 2 (7V + l) + ( 5 li2)1 (7V + l) 



|c 2 (5 2 (iV + 1) + S 2 (N + 1)) + (^(N + 1) 



1 



2S 1>3 (N + 1) + S 3;1 (N + 1) + S 2 , 2 (N + 1) 



N + l 

-C 2 5 2 (iV+l)-2C 3 5 1 (iV+l) 



(B.45) 



(B.46) 



27 



By summation one obtains from (B.43 
/ ln(l — x)Li 2 (l — x)\ 



-B.46) 



M 



x — 1 



M 



M 



'm 2 (l - x)U 2 (l - x) 
v x — 1 

^m(l - x)Li 3 (l - x) 



x - 1 



M 



'ln(l - x)S li2 (l - x) 



x- 1 



;-i) jv m 



(-l^M 



ln(l - x)Li 2 (l - x) 
x+ 1 



ln 2 (l -x)Li 2 (l -x) 
x+ 1 



(N) = 2S 1 , lt2 (N)+S 1 , 2 ,i(N)-2<; 2 S 1 , 1 (N) 

= -S 1 (N)S 2 , 1 (N)-S 3il (N)-S 2 , 2 (N) 
+2S 3 (N)S 1 (N) + 2S 4 (N) 
+2(S 2 (N)-{ 2 )S 1 , 1 (N) 



(N) 



(TV) 



----- -6S , i 1 i,i l2 (JV) - 4S , i,i )2l i(iV) - 25 1)2il)1 (7V) 

+6C 2 5 , i,i 1 i(iV)+2C35 , i 1 i(N) 
-3S , i 1 i,i, 2 (JV) - 5 lilj2il (iV) + 3C 2 S'i,i,i(iV) 



(TV) 



(AT) 



(-l^M 



ln(l — x)Li 3 (l — x) 
x+ 1 



(N) 



'-1) N M 



ln(l-x)S 1>2 (l-x) 
x + 1 



(AO 



(B.47) 



(B.48) 



-C 3 Si,i(A0 
= -S 1 (N)S 1X2 (N) + S 1 ^ 2 (N) 

S 2 (N)S li2 (N) - S 3 , 2 (N) - S 1A (N) - Si,i, 3 (A0 
+3C 2 Si ) i 1 i(W) - Si,i(A% (B.49) 

= 2S hl:3 (N) + 5 1)3 ,i(AT) + S\ :2j2 (N) - ( 2 S h2 (N) 

-2( 3 S hl (N) (B.50) 
= 2S_ 1;1 , 2 (A0 + S'_i )2) i(JV) - 2C 2 S- ltl (N) 

1\ 3 



-^ + ^n 4 (2) + 3Li 4v2 



ln 2 (2) 



(B.51) 



(iV) = -6S , _i,i i i, 2 (iV) - 4S-_ lil)2il (JV) - 2S-_ 1)2)1)1 (7V) 
+6C 2 ^-i,i,i(A r ) + 2C 3 5_ 1 , 1 (iV) 
-^C 2 Cs - ^2 ln 3 (2) + gc 2 2 ln(2) - 1 ln 5 (2) 



4 



20 



10 



(B.52) 



-3S'_ 1)lil)2 (JV) - S- ltl , 2>1 (N) + 3( 2 S^ 1X1 (N) 

-C 3 ^ 1 , 1 (iV) + -C 2 C3 + ^C 2 ln 3 (2) 

_^ C 2 ln(2) l c ln2(2) _ 29 + 1 ln 5 (2) 



-4Li fi 



(B.53) 



2S_ 1;1 , 3 (A0 + S^ 3A (N) + S^ 2>2 (N) - ( 2 S_ 1>2 (N) 
-2{ 3 S- 1A (N) 

-f C2C3 + ln 3 (2) - ^Cl ln(2) + ^3 ln 2 (2) 
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+ 



159 



64 



(B.54) 



The above expressions depend on harmonic sums, for which we derive the following relations. 
They are used to derive the asymptotic representations only. 



S-i,i,i(N) 
Si,i, 2 (A0 



5-1,3,1 (AO 



S 1 (N)S- 1 , 1 (N) + 5_ 2;1 (7V) + S- li2 (N) 
-S 1 (N)S li _ 1 (N) - S 2 ^(N) - S h _ 2 (N) 



[-1) n \f 3 (N) + ln(2) [S^N) - S ltl (N)] + \ ln 2 (2) \S X {N) - S^N)} 
1 



6 



In 3 (2) 



= -F 7 (N) + C 2 S , i,i(JV) - Cs-S^JV) + -C 2 5 



(B.55) 

(B.56) 
(B.57) 



:-l) N F 8 (N) + ( 2 S- 1A (N) - CsS^(N) - Li 4 Q + Acl - ^3 ln(2) 
-iC2ln 2 (2)-lln 4 (2) 



(B.58) 



(-l^M 



251,2 ( -x) + ln(l + x)Li 2 (-x) 



1 + x 
iC3-^C 2 ln(2) 



(iV)-ln(2) [S 1:2 (N) - S lt _ 2 (N)] 



-2^1,-1^) + 
-^C 3 ln(2) + C 2 ln 2 (2)--ln 4 (2) 

-F 12 (N) + C 3 Si,i(iV) - ^i(iV) + 2C 5 - C 2 Cs 
1 



[S^N) - S^(N)} + ^C! ~ 3Li 4 (i 



(B.59) 
(B.60) 



-F 13 (iV) - 2S- lil)3 (iV) + ( 2 S lj2 (N) + 2{ 3 S ltl (N) - fS^N) + ( 2 ( 3 - 



S-i,i,s(N) = (~1) N F 12 (N) + C 3 S-i,i(iV) - $:S-t(N) + 2 



Li,(i)+LU Q) ■ 



+^ln 5 (2) - ic 2 m 3 (2) + ^(2) + -^Ca - ^C 5 



(B.61) 



(B.62) 



= S^iN) + S 1 (N)S^ 3 (N) + S- 1A (N) - S-i(N)Si t3 (N) - S lj ^(N) (B.63) 



= (-1) 



iV+l 



M 



ln(l + x)Li 3 (-:r) + U\{-x)/2 
1 +x 



(N) + \n(2) [S , i,_ 3 (iV) - S 1;3 (A0] 



+ (4 C ' + I C3hl ( 2 )) " " ^ 2 Sl,- 2 (iV) + ^3^1,-1 W 

1 1 1 ~\2 e ) 1 

+-C|ln(2) + -C 3 ln 2 (2) - -C 2 ln 3 (2) - —(5 + -hi 5 (2) 



+2 



Li,(^)+LU Q) 



— T^C 2 C3 

lb 



5-l,2,2(N) = 



(-1)^14^) - 25_ 1)1)3 (iV) + C 2 5_i )2 (iV) + 2C 3 5_ lil (7V) - |5-i(A^) 
_4Li5 G) + ^ ln5(2) " <^ 2 ln3(2) " S C2 ln(2) + C2Cs + Y Cs 



(B.64) 



(B.65) 
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S-i,i,i, 2 (AO = -(-lTF 18 (N) + ( 2 S- 1 , lil (N)-( 3 S- lil (N) + -(tS- 1 (N) 



2Li 5 - -ln(2)Li 4 



-lln 5 (2) + ic 2 ln 3 (2)-^C 2 2 ln(2) 



7 1 



(B.66) 



7V+1 



M 



3S 1>3 (x) + 2S 1 , 2 (x) ln(l - g) + ln 2 (l - x)Li 2 (g)/2 

1 + :r 



(iV) 



+C 2 S-i,i,i(AO " 4C 2 C,3 - ^C 2 ln 3 (2) - ^C 2 2 ln(2) + ^C 3 ln 2 (2) - ^5 



16 



12 



40 



16 



64 



13 



5-1,2,1,1 (AO 



120 
(-1) W M 



ln 5 (2) + 31n(2)Li 4 Q + 2Li 5 Q 
ln(l-x)S' 1)2 (a;) + 3S , li 3(x 



(B.67) 



1 + x 



(N) + ( 3 S_ 1A (N) 
11 



+ y C.C3 + ^C 2 ln 3 (2) - ^C 2 2 ln(2) - ^Cs ln 2 (2) - ^ ln 5 (2) 



+ g C5 -31n(2)Li 4 (I)-4Li 5 Q 

The above relations require the function 

~ln 3 (l -x ~ 



M 



1 +x 



(B.68) 



(B.69) 



which is related to M 

Si,i,i,-i(A0 = ' 
S-i,i,i,i(N) = 



ln 3 (l +x)/(l +x) (AT) by 



-1) W+1 -M 
; 6 

-lf +1 -M 



ln 3 (l - x) 
1 +x 

ln 3 (l + :r) 
1 + x 



(iV) - Li 4 



(iV) + ln(2) [Si,i,_i(JV) - Si,i,i(JV)] 



5_ 1 , 1 , 1 (7V) + 5 1)1 ,_ 1 (7V)--5 1 ,_ 1 , 1 (7V) 



(B.70) 
(B.71) 



~ [^ 2 ,-i,i(iV) + Si,- 2 ,i(A0 + S , i_i l2 (iV)] + ^[^_ 1)li2 (7V) + 5_ 2)lil (7V) 



+5 , -l,2,l(A 7 ) + Sl,l,-2( N ) + 5 2,i,-i(A0 + S lt2 ,-i( N ) 
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The sums associated to the index sets {1,1, —1} and {1, 2, —1} are represented by F 2 (N), resp. 
Si,2,-i(N) and S-± t i i2 (N), using the algebraic relations [15]. We use the subsidiary functions 



#i(A0 = 


M 


~2S lt2 (-x) +ln(l + x)Li 2 (-x)" 


1 + x 


H 2 (N) = 


M 


"ln(l + x)Li 3 (-a;) + Li|(-a;)/2 


1 +x 


H 3 (N) = 


M 


1 + x V ; 



(N) 
(N) 



(B.73) 
(B.74) 
(B.75) 



to obtain the asymptotic representations, cf. Sect. B.4. The other sums are two-fold and their 
explicit representations were given in [2] and in the previous sections. 
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Changing the arguments from z to 1 — z in the numerator functions leads to relations like 
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M 



M 



M 
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z±l 
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(TV) 



(TV) 
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-M 
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(TV) 



(TV) 
(TV) 



(B.81) 



which are used to express the asymptotic representations, Sect. B.4. 
B.3 Recursion Relations 

The following recursion relations hold for the basic functions : 

i>(i + z) = 4>{z) + - 

z 

^ n \i + z ) = ^ n \z) + {-iy 
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(B.82) 
(B.83) 
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8(2 + l) 2 2(z + l) 3 



1 



d 



r ln 2 (2) 



2(* + 1) dz 



z+1 z+1 



FAz + l) 



(B.107) 



cf. (7.3). 

The recursion relations can be obtained from the integral-representations of the Mellin 
transforms. Any point in the analytic region of the basic functions is connected by the recursions 
which have to be applied until the asymptotic region \z\ > £ asympt with -2 asympt ~ 15. ..20 is 
reached. 



B.4 Asymptotic Representations 

The asymptotic representations of functions contributing to the single harmonic sums and the 
basic functions Fi(z) and Fi(z) are given in the following. We maintain contributions up to 
0(1/ z 19 ), to reach double precision accuracy in a fully analytic representation. Higher order 
terms can calculated if needed in more special numerical applications. 

1 00 R 
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^ 2 \z) 



n-l)\ n\ ~ D (2k + n 



z n 2z n+1 ^ ' K (2k-1)\ Z 2k+n 



(B.109) 



,,,1111 11 11 11 11 691 

xjj(z) ~ \n(z)---- — — + — — - — — + — — - — — + 



2 z 12 z 2 120 z 4 252 z 6 240 z 8 132 z 10 32760 z 12 
1 1 3617 1 43867 1 174611 1 n f 1 A 
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The asymptotic representations of the subsidiary functions #1,(2,3) (z) are given by : 
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In the above Tk denotes the tangent numbers or Euler zigzag numbers, [27,60], with the gener- 
ating function 



oxp(,r) -rxp( ,r) ^ (-l) fc 1 2fc _ 1 



E 



exp(x) + exp(— x) (2 A; — 1)! 

They are related to the Bernoulli numbers by 
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C Some Integrals of w = 5 

In this appendix a series of integrals used in the present calculation is presented. Other integrals 
of this type were given in [63]. 
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